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Abstract. In this paper we consider a large class of many-variable polyno- 
mials which contains generalisations of the classical Hermite, Laguerre, Jacobi 
and Bessel polynomials as special cases, and which occur as the polynomial 
part in the eigenfunctions of Calogero-Sutherland type operators and their 
deformations recently found and studied by Chalykh, Feigin, Sergeev, and 
Veselov. We present a unified and explicit construction of all these polynomi- 
als. 



1. Introduction 

In this paper we discuss the construction of symmetric polynomials which arise 
as eigenfunctions of exactly solvable quantum many-body systems of Calogero- 
Sutherland type [6,51]. In particular, we demonstrate that a particular construc- 
tion method found and studied by us in simple special cases [20,30] can be naturally 
extended to the full class of such model. Our results provide a unified approach to 
many-variable generalisations of the classical Hermite, Lagueree, Jacobi and Bessel 
polynomials. We also show that it is natural to generalise our approach to a larger 
family of polynomials related to a deformation of the Calogero-Sutherland type 
systems found and studied by Chalykh, Feigin, Sergeev, and Veselov [7,46,47]. 
This allows us to derive, for each such polynomial, a family of different representa- 
tions labeled by a pair of non-negative integers (M, M). Moreover, one can choose 
(M, M) such that the representation is simplest, i.e., a linear combination of certain 
explicitly given polynomials with the least number of terms. These results are non- 
trivial already for the non-deformed polynomials. To mention a specific example, 
we find that a Jack polynomial [49] has a (M, M) representation if there exists a 
pair of non-negative integers (m, m) such that the Young diagram corresponding 
to this Jack polynomial can be covered by the union of the two rectangular Young 
diagrams of size M x m (i.e., M rows and m columns) and rh x M, respectively. 
Moreover, the simplest such representation is obtained by minimizing M + M. For 
example, the simplest possible representation of the Jack polynomial corresponding 
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to the Young diagrams 



» 



and 



is for (M,M) = (1,0), (0,1), and (2,3), respectively, and in the first two cases 
our representation actually consists only of a single term. The main tool in our 
construction are certain identities which, to each (deformed) Calogero-Sutherland 
system, relates a system of the same type. In special cases, these identities are well 
known and related to a duality of the Jack polynomials (see, e.g., [39]), or, more 
generally, the algebra homomorphism which maps Jack polynomials to so called 
super Jack polynomials (see, e.g., [48]). We deduce these identities in a unified 
manner, and thus unify as well as extend, previously known results. 

While we mention these examples in the beginning to give an impression of the 
kind of results we obtain, the emphasis in our presentation is on our construction 
method. To avoid overloading our formulae we explain this method in detail for 
the non-deformed polynomials and the representations (M, 0). The subsequent 
extension to the general case is then remarkably easy. 

1.1. Calogero-Sutherland type models. As mentioned above, the polynomials 
studied in this paper play an important role in physics since they occur in promi- 
nent exactly solvable quantum mechanical models defined by partial differential 
operators of Calogero-Sutherland type. These latter operators play a central role 
in our construction, and we therefore proceed to explain their relation to the clas- 
sical Hermite, Laguerre, Jacobi and Bessel polynomials. As elaborated below, to 
each such sequence of one- variable polynomials {p n : n = 0, 1, . . .} there exists a 
quantum mechanical model defined by a Schrodinger operator 

h = -dl 



(1) 



V(x), 



with a particular potential function V such that its eigenfunctions are of the form 



(2) 



1p n {x) = 1p (x)Pn(z) With Z = z(x) 



for particular C 2 functions ipo(x) and z ( x ) (% is a real variable, and we use the 
notation d x :— Jj). Many- variable generalisation of these polynomials appear in 
the following generalisation to an arbitrary number N of particles: each one-particle 
Schrodinger operator |T]) allows for a many-body generalisation 



(3) 



H 



N 



N 

E 

3=1 



j<k 



with W a particular two-body interaction potential such that iJ/v has eigenfunctions 
\l/ n which are labelled by integer vectors n = (n\, . . . , ri/v) and are of the form 

(4) * n (x) = * (x)P n (z) with Zj = z( Xj ), 



N 



(5) 



* (x)=n^o(x j )n(^-^r 
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where the functions P n are certain symmetric polynomials in TV variables z = 
(zi, . . . ,Zn) generalising the polynomials p n . To avoid certain technicalities and 
since it includes the cases of main physical interests we will assume k > 0, even 
though many of our results actually hold true also for other values of k. We mention 
already at this point that although the symmetry of the polynomials P n implies 
that a complete set of eigenfunctions can be labelled by partitions, i.e., sequences 
n = A = (Ai, . . . , Aj\r) of non- negative integers Aj such that Ai > ■ ■ ■ > Ajv, it is in 
our approach not always natural to make this restriction, as will become evident 
below. These operators Hn define exactly solvable quantum-many body systems of 
Calogero-Sutherland type [6,22,42,50]; see Table Q] for a list of well-known examples 
(note that we have added constants to the potentials V to simplify certain formulae 
later on). 

1.2. Deformed Calogero-Sutherland operators. A further interesting gener- 
alisation of the one-particle Schrodinger operators |T]) is the following class of dif- 
ferential operators in two sets of variables x = . . . , xn) and x = (x%, . . . , Xfj), 
where N and N are arbitrary non-negative integers: 

N N 

(6) H N f! = ]T {-d 2 x . + V{ Xj )) - J2 *(-dl + V"(a:j)) 

j=l ,7=1 

+ K(K-l)YlW(Xj,Xk) + V W{Xj,X k ) 

i — < K t — < 

j<k J<K 

+ (1- K)^2W(Xj,X K ) 

with a potential function V of the same kind as V but with different parameters, as 
specified in Table [TJ These differential operators are natural generalisations of the 
Schrodinger operator ([3]) in that they also have polynomial eigenfunctions with cor- 
responding eigenvalues that can be computed explicitly; see Section [4] for a precise 
formulation of this result. In this paper we present an explicit construction of all 
these many- variable generalisations of the above mentioned classical polynomials. 
We obtain the results by a unified treatment making no reference to special cases, 
as explained below. 

1.3. Construction method. To explain the nature of our results we now dis- 
cuss the standard Calogero-Sutherland cases N = in more detail. Our approach 
is based on a set of remarkable identities which provide us with particular sets 
of symmetric polynomials on which the action of the differential operators ([3]) is 
particularly simple. These latter polynomials can be defined by the following ex- 
pansion: 

(7) Ul ^ k -f " w ;' w f = E /n^r'< a • • - v, 

valid for \wn\ > |iUjv-l| > ••• > \wi\ > niaxfe(|^fe|)- We mention that these 
polynomials were first defined in [30] by a certain contour integral. However, using 
Cauchy's theorem it is easy to verify that the two definitions are equivalent. We 
explicitly construct series representations for the polynomials P n in which are of 
the form P n (z) = J2 m u n(ni)/ n (z) and, in addition, we obtain certain completeness 
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results; see Section 3 for the precise statements. We will also derive a partial 
generalisation of these results to the differential operators H N see Section 4. We 
thus generalise, unify, and extend various results which were known before only in 
special cases [20,31]. 

We now explain in which sense our approach is unified and give a precise def- 
inition of the classical polynomials we consider: each such sequence consists of 
polynomials p n of order n = 0, 1, . . . such that they are common eigcnfunctions of 
a differential operator 



(8) 
with 



-ij- 1 h^ = a(z)d 2 z +P(z)d z 



(9) 



a{z) = a-iz 1 + a.\z + ao and (3{z) = (3\z + /3q 



for some real coefficients ctj and f3j ; see Table [2] for the polynomials corresponding 
to our examples in Table [TJ It is interesting to note that the Schrodinger operator 
h in (fT]) with the eigenfunctions ([2|) can be fully characterised by these polynomials 
a and (3: it is straightforward to verify that 

(2/3(z) - a'(z))(2(3(z) ~ 3a'(z)) 



(10) V{x) = v(z(x)), v(z) 



16a(z) 



±a»[z) + \(f{z) 



V{x) 



W(x,y) 



V(x) 



I 

II 
III 

IV 

V 
VI 

VII 



2 2 
UJ X 



cosh 2 (:r) U 



. ,2„2 , a(a-l) 
tj-> X 1 

-w(l + 2a) 



n(a-l) _2 
— 5 — a 



a(a-l) . b(fr-l) 
4sh?(§) ' 4cos 2 (§) 
(n+b) 2 



W 2 e~ 2a; 



-w(l + 2c)e- ;r + c 2 



2 

(z-'S) 2 
1 

2sin 2 (^) 
1 i 



2sinh 2 (^) 2cosh 2 (^i) 
2 , 2 



2sin 2 (^i) 2sin 2 (^i«) 



+ 1 



2s i„ 2 (2L=Ji) ' 2sin 2 (£S) 



2sinh 2 (^i) 



(K 2 - 1)/k 4 
__c(c+iL 

cosh 2 (2; ) 



w z x z H 

-w(l + 25) 



q(a-l) ~2 
— — a 



q(q-l) , b(b-l) 
^C£1 "T" J777727 



4sin 2 (f ) 4cos 2 (f ) 
(5+b) 2 



c^e" 2 * 



-0(1 + 2c)e-* + c 2 



Table 1. Examples of external- and two-body interaction poten- 
tials such that the differential operators ([3]) and ((6]) have poly- 
nomial eigenfunctions. The parameters in the last column are: 
Cj = —oj/k, c — —(2c + k + 1)/(2k), a = —(2a — k — 1)/(2k) and 
b = -{2b - k - 1)/ {2k). 
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(the prime here and in the following indicates differentiation) with z(x) a solution 
of the differential equation 

(11) z'(x) 2 = a(z(x)) 
and 

(12) i> (x) = e - w W x » 
with w(z) any solution of the differential equation 

y ' Aa{z) 

Moreover, the potential functions W, ensuring that the differential operators Hn 
in §3§ and H N ^ in ^ have polynomial eigenfunctions, are given by 

a(zi) + OiiZi) 

(13) W(xi,X2) = — -, tt. a 2 with Zj = zixj). 

{zi- z 2 ) 2 

The constant a 2 is subtracted to simplify formulae later on, and V is given by a 
formula as in ifTU)) but with (3{z) replaced by 

0(z) = [(1 + k)c/(z)-0(z)]/k. 

It is interesting to note that the exact eigenvalues of the differential operators Hjsr in 
([3]) and H N # in ([6]) are determined by the leading coefficients of these polynomials 
a{z) and @(z). For example, the exact eigenvalues of Hn corresponding to the 
eigenfunctions (|3|) are 

N 

(14) E a = E Q - (^njirij - 1) + (2Ka 2 {N - j) + (3i)n 3 ) 

j'=i 

with 

(15) E = -^N(N 2 1) - + N{N 1} . 

This result, well-known in many special cases (see, e.g., [14]), will naturally emerge 
from our construction in Section^ We are now in a position to state more precisely 
in which sense our approach is unified: in our construction information concerning 
the different cases will enter only via the two polynomials a and (3. We also stress 
that the special cases listed in Table [2] are only intended as examples, and that our 
approach is for arbitrary polynomials a and (3 of the form ([9|) . 

1.4. Special cases. To put this general scheme into perspective we now discuss 
its well-known special cases summarised in Tables [T] and [5] The case a = 1 and 
j3(z) = —2luz corresponds to the exactly solvable many-body generalisation of the 
quantum harmonic oscillator introduced by Calogero [6] (case I in Tables [I] and 
[5]), while the Sutherland model [50,51], generalising the quantum model of a free 
particle on a circle, is obtained by setting a(z) = — z 2 and (3(z) — ~z (case II). 
These models correspond to the Hermite polynomials and the ordinary monomials 
Pn(z) = z n , respectively. The many-body models associated with the Legendre, 
Gegenbauer and Jacobi polynomials (cases IV, V and VI) are identical with 01- 
shanetsky and Perelomov's B^- and BCAr-variants of the Calogero- and Suther- 
land models; see [42] and references therein. An exact solution to the many-body 
generalisation of the Morse potential (case VII), corresponding to the generalised 
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ipo (x) 


z(x) 


Pn(z) 


a(z) 


P(z) 


I 


g-wa; /2 


a; 


H n («/ujz) 
(Hermite) 


1 


—2tuz 


II 


1 




z n 


-z 2 


—z 


TTT 
111 


cosn [x) 


i sinh(x) 


C n (Zj 

(Gegenbauer) 


1 i -,2 
— 1 + Z 


(1 — zcjz 


IV 




a; 


L n (ujz) 
(Laguerre) 


4z 


z + 4a — 4cjz 


V 


sin a (x) 


cos(x) 


(Gegenbauer) 


1-z 2 


-(1 + 2a)z 


VI 


sin a (f)cos b (§) 


cos(x) 


(Jacobi) 


1-z 2 


6 — a— 
(l + a + b)z 


VII 


exp (~u;e~ x — ex) 


e x 


y n (z, 1 - 2c, 2uj) 
(gen. Bessel) 


z 2 


2u + (1 - 2c)z 



Table 2. Exact solutions of the one-body Schrodinger equation 
associated with classical polynomials, as described in (JTJ) and @. 
Given are also the associated polynomials a and (3; see (|8T)- (fT2|) . 
More details about the polynomials p n can be found in [1] (cases I 
and III- VI) and in [16] (case VII). 



Bessel polynomials [26], was first given by Inozemtsev and Meshcheryakov [22]. A 
further study of this case can be found in [19]. There exist various other interesting 
special cases which we have not explored in detail, but they can all be transformed 
to the cases mentioned above by a rescaling and a translation of the variable z. 
The deformations of the Calogero-Sutherland models defined by H N ^ in (J6j) were 
found and explored by Chalykh, Feigin, Sergeev, and Veselov; see, e.g., [7,46-48]. 
We will refer to these differential operators H N ^ as deformed Calogero-Sutherland 
operators. 

1.5. Related previous work. We now discuss the relation of our results to pre- 
vious work in the literature. The many-variable polynomials corresponding to the 
eigenfunctions of the Calogero-Sutherland type models in cases I- VI in Table Q] 
have been extensively studied in the mathematics literature; see, e.g., [11,39] and 
references therein. We mention, in particular, Heckman and Opdam's root sys- 
tem generalisation of the Jacobi polynomials [21] and the work of Baker and For- 
rester [3], van Diejen [9], and Lassalle [33-35], as well as of Macdonald [38], on 
many-variable generalisations of the classical orthogonal Hermite, Laguerre and 
Jacobi polynomials. A particularly well-studied case are the so-called Jack poly- 
nomials corresponding to the Sutherland model (case II). Explicit formulae for 
the Jack polynomials were recently obtained by Lassalle and Schlosser [37] (see 
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also [36]) by inverting a so-called Pieri formula. For very particular values of the 
integer vector n = (ni, . . . , tin) or a low number of variables explicit series expan- 
sions of the Jack polynomials were obtained by Stanley [49]. In addition, there exist 
expansions of the Jack, as well as certain other many-variable classical orthogonal 
polynomials, which are of a combinatorial nature [8,10,27,39]. These combinato- 
rial results also include a formula resulting from the use of Sutherland's original 
solution algorithm [51]. We also mention operator solutions of the Calogero- and 
Sutherland models obtained in [5,24,25,52,53] and [32], respectively, as well as in- 
tegral representations of the Jack polynomials; see [2,40,41] and references therein. 
Our results in this paper provide an explicit construction of the many-variable 
polynomials P n . These results were recently announced in [17], and they gener- 
alise those previously obtained by us for the Sutherland- [30] and the Calogero 
models [20] . We mention that our unified treatment of all the cases in Table [1] is 
different from the one based on root systems (see, e.g., [42]), and that it has been 
previously used by Gomez-Ullate et.al. [14] to obtain the energy eigenvalues in (fM]) 
and (|15p by purely algebraic means. We also mention that our point of view is 
closely related to the theory of quasi-exactly solvable Schrodinger operators; see, 
e.g., [4, 15]. Moreover, the identities which are the key to our results are stated in 
Corollaries 12.21 and 12.31 below, and they were known before only in special cases: 
an important special case of these identities for the Sutherland model (case II) is a 
consequence of a well-known result on Jack polynomials which, to our knowledge, is 
due to Stanley (see Proposition 2.1 in [49]), and a generalisation of the latter to the 
deformed case and other non-deformed cases can be found in [48] and [13,39,45], 
respectively. These identities relate Schrodinger- or deformed Calogero-Sutherland 
operators with different parameter values and are natural quantum analogs of the 
Backlund transformations for the classical Calogero-Moser models first found by 
Woijekowski [54], as discussed by Kutznetsov and Sklyanin [29]. We obtain these 
identities as specialisations of a particular identity stated in Proposition 12 . 1 1 which 
has the natural physical interpretation of giving the groundstate of a generalisation 
of the Schrodinger operator where the particles are allowed to have different 
masses rrij and with specific mass dependent external potentials V m . (xj ) . This is 
a powerful result of independent interest, and to our knowledge it was previously 
known only in the special cases I, II and IV in Table [T] [12,20,44]. We finally 
mention an integral representation of the Jack polynomials recently obtained by 
Kuznetsov et.al. [28] using a separation-of-variables approach which is also based 
on the identity in Corollarv l2.2l 

1.6. Plan of the paper. In Section[2]we derive and discuss the identities which are 
the key to our solution method. Section[3]contains our explicit construction of series 
representation of all many-variable polynomials P n determining the eigenfunctions 
(Q of the Schrodinger operators ([3]). In Section [4] we extend our construction of 
polynomial eigenfunctions to the deformed Calogero-Sutherland operators ((6]). We 
conclude with a few remarks in Section [5] 

1.7. Notation. We denote by Z, N, No and M the sets of all integers, positive 
integers, non-negative integers and real numbers, respectively. We shall say that 
an integer vector n = (n%, ri2, . . .) is of length N, denoted £ (n) = N, if N is the 
smallest non-negative integer such that the parts rij of n are zero for j > N. The 
symbols n and m will in most instances denote integer vectors, except in Section [2] 
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where the symbols m,j are used for real 'mass parameters' rather than for parts 
of an integer vector m. We will use the symbols A, /J, and u to emphasise that a 
particular integer vector is a partition. We will also use the notation 

Y s _ si s N 

where x = (xi, . . . , xn) and s = (si, . . . , sat). Moreover, for m 6 7L M and n 6 Z N 
we will write (m, n) short for (mi, . . . , mju, n\, . . . , tin). 

2. Identities 

In this section we present and prove a particular identity associated with the 
one-particle Schrodinger operator ([1]). The identities underlying our construction 
of eigenfunctions of the differential operators in ([3]) and ([6]) are stated in Section l2~2l 
and obtained as special cases of this more general results. Throughout the section, 
we assume the polynomials a and j3 in ([9]) fixed, and z(x) and W{x\, x%) are as in 
£3} and 

2.1. Source identity. The following identity can be interpreted as providing the 
exact groundstate of a generalisation of the Schrodinger operator (|3|) where the 
particles can have different masses rtij. However, we will allow these parameters 
rrij also to be negative, and we will use this identity as a source from which we obtain 
various other identities as special cases. To state this result we find it convenient 
to denote this generalised Schrodinger operators as H, the particle coordinates as 
Xj , and the particle number as A/". 

Proposition 2.1. Let 

M 

(16) H = — (-9% + V mj (X,)) + ^Km* - l)(m, + m^WiX^Xk) 

j=i 1713 j<k 

with Tn j arbitrary real and non-zero parameters and 

(17) V m (X) =v m (z(X)) 
with 

_ (2/3 m (z) - a'(z))(2p m (z) - 3a' (z)) 1 I,, 
v m {z)- (z)+ 2 p m (z), 

where the prime indicates differentiation with respect to the argument z and 

(18) (3 m (z) = m(3{z) + m)(l - nm)a'{z). 
Furthermore, let 

AT 

$ (x) = n (x,) H(z k Zj r m ^ 

j=l j<k 

with Zj — z(Xj) and 

(19) ipo, m (X) = e -™»(«W) 
with w m {z) any solution of the differential equation 

a'(z)-2/3 m {z) 



4a(z) 
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(H - So) $o = 



Then 
(20) 

with the constant 

(2D £o = ~^(H 3 - l- 3 D - "fo-t^W - I- 2 ! 



where 
(22) 



AT 



l m "l : =E m " »= 1,2, 3. 



The m-dependence of the external potentials and one-particle groundstate eigen- 
functions for our examples is given in Table[3]with m-dependent parameters dehned 
in the table captions. 

Remark 2.1. Note that (3\ and (3q always refers to the coefficients of the polynomial 
(3(z) as dehned in @, and are not to be confused with f3 m {z) for m = 1 and m = 0. 
There should be no danger of confusion since we always write (3(z) for m (z) if 
m = 1 and assume m^O. 



V m (x) 



I 

II 
III 
IV 

V 
VI 
VII 



jK 2 m 2 (m — l) 2 



c m (c ? + l) 
cosh 2 (a;) 

, ,2 2 , a m (am-l) 
W m X H -3 



a rn {a rn — l) 
sin 2 x 

a m (a m -l) j_ b m (b m -l) 
4 sin 



■ u m (l + 2a m ) 



^(f) - + 4cos 2 (f ) 4V a mTO m J 

2x - w m (l + 2c m )e- a: + (? m 



e -u, m x 2 /2 

cosh _Cm (x) 

p — u m x 2 /2 rJ ,a rn 



sin am (x) 
sin a '"(f)cos b '"(§; 



Table 3. 'Mass' dependence of the external potentials and one- 
particle groundstates in our examples, according to Proposi- 
tion 12.11 The m-dependence of the parameters is as follows, 
u> m = mui, c m — mc — ^Km(m — 1), a m = ma + ^Km(m — 1), 
and b m — vnb + ^Km(m — 1). 



Proof of Proposition \2.1\ We will show by straightforward computations that 

N 1 

(23) n = 12— QjQJ + £o, 



3 = 1 



where 



Qf=Tdx j +V j , V i (x) = 
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Since Q } - $o = for all j, this will prove the identity in (j20|) . 
To prove (|23[) we compute 

Vj (x) = -w' mj {z{X 3 ))z' (X, ) + J2 Krn jmk : ' ' A ' ' 



and thus 



E = - E + wi + w 2 + w 3 

j — 1 i — l 3 



with 



M 1 / 1 



J=l 3 



the one-body terms, 

2 = SS""" 1 <K """" " '(^^ — ^ — . 

the two-body terms, and 

a(Zj) 



W3 = E E K 2 mjm k m e - 



(Zj - Z k )(Z, - ZA 

the three-body terms; we inserted the relations z'(Xj) 2 = a(Zj) and z"(Xj) = 
ha'(Zj) and used the short hand notation z(Xj) = Zj. 

Symmetrising the sum defining W3 and inserting a(Zj) = ^2p=oCtpZj we can 
write 

^ Z]{Z k - Z t ) - Zj{Z 3 - ZJ) + ZjjZj - Z k ) 
W 3 = 2^ 2k m 3 m k m e }^a p — — — — — — , 

and observing that the fraction on the r.h.s. is identical with and 1 for p = 0, 1 
and p — 2, respectively, we find that W3 is equal to the constant 

W3 = 2K 2 mjm k mia.2- 

j<k<£ 

Symmetrising the sum defining W2 and inserting ^a'(Zj) — 2w' m . (Zj)a(Zj) = 
(3 m AZj), following from the second equation in (|19|) . we obtain 

... \ - (< m k a(Zj) + mja(Z k ) m k /3 mj (Zj) - mj/3 mk (Z k )\ 

m = § K { {Kmjmk ~ 1} &-z k y + z—z k J ■ 

We now decompose W 2 into two parts as follows: 

m = W 2 ,o + W 2 ,i, 

where 

K , ,i, x a (Zj) + a(Z k ) 

W 2 ,o = 2^ ~( Km j m k ~ l)(n»j + m k )- 



2 1 (Z,--Z k ) 2 
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and 



j<k \ 



\— | i ,a{Zj) — a{Z k ) 
yv 2 ,i = y j k\ -{Kinjirik - l)(m k - rrij) — 



, m k /3 m] (Zj) - mj(3 mk (Z k ) 



Recalling Q13p we find that the symmetric part gives us the two-body terms we 
want, up to a constant, 

W 2 ,o = ^2 -z( Km i m k ~ + TO fc ) (W(xj,x k ) + a 2 ) ■ 

j<k 

Inserting (fT8|) and © a straightforward but somewhat tedious computation shows 
that the terms in the antisymmetric term add up to a constant, 

2~ 



W 2 ,i = ^2 K [—(KmjTrik + l)(rrij + m k ) + (/3i - a 2 (l + K))mjm k 
j<k 



Inserting the second equation in (| 19(1 and (|18p a simple computation shows that 
the one-body terms are identical with 



Wi = z2— v ™A x j) 



with V m defined in (fl7| . Collecting all terms we obtain the identity in (|23|) with 
the constant 

£ = — 2k 2 a 2 rrijm k mi — K 2 a 2 ^2 m j m k{ m j + m fc) 

j<k<l j<k 

- k(Pi - a 2 (l + k)) ^ mjin k . 

Using the notation in (|22| we find by straightforward computations that this con- 
stant is identical with the one given in (|2"T|) . □ 

The physical interpretation of TL in (|16p as Schrodinger operator of a quantum 
many-body system requires that it defines a self-adjoint Hilbert space operator 
bounded from below. This is the case under certain obvious restrictions on param- 
eters. We now discuss this Hilbert space structure for the different cases listed in 
Table [TJ but our discussion will be brief since this aspect does not play any role for 
our construction in the following sections. 

In cases I, III and VII the relevant Hilbert space is L 2 (]R A ^), in case IV we have 
instead L 2 (K+), while £ 2 ([-tt, tt] m ) is associated to case II, and i 2 ([0, tt]^) to 
cases V and VI (the weight function is in all these cases constant and equal to 
1). It is obvious in all these cases that in the proof above is the Hilbert space 
adjoint of QJ (on suitable domains), and thus, if all rrij are positive, that Ti. in (|23|) 
defines a unique self-adjoint operator via the Friedrichs extension (see, e.g., [43]) 
with <!>o as groundstate provided that $o is square integrable. In particular, the 
Hn in ((3]) define, under obvious restrictions on parameters, self-adjoint operators 
bounded from below. However, this is not the case for the deformed operators ((6]). 
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2.2. Important special cases. We proceed to discuss the special cases of Propo- 
sition 12.11 which underlies our construction of eigenfunctions of the differential op- 
erators © and © in Sections [3] and [U respectively. Setting Af = N, rrij = 1 and 
Xj = Xj for j = 1, . . . , N we obtain as a first special case the following: 

Corollary 2.1. The function ^o, as defined by §E§ and (| 1 2(1 . is an eigenfunction 
of the Schrodinger operator ([3]) with corresponding eigenvalue Eq ■ 

The remarkable identity underlying our solution method in Section [3] is obtained 
from Proposition 12.11 by choosing J\f — 2N, nij = 1, Xj = Xj, m^ + j = — 1, and 
Xn+j = Vj for j = 1 , . . . , N. We observe that Ti in (fl6|) then becomes a difference 
of two Schrodinger operators (J3j> . Denoting $o as F and £q as Cn we obtain the 
following: 

Corollary 2.2. Let -ffjv(x) = Hm be the operator in J3J), 

N 

(24) H [ N \y) =E(-<+ + K ( K - X ) E ^(%. Vk) 



■j=i j<k 



and 



F(x, y ) = y_ n fc K-^) w 

J— 1 J' 

wii/i = z(xj), Wj — z(yj), and V—i(x) and -0o._i(x) as in (j 1 T[) and (|19|) /or 
m = — 1. T/ien 

(25) (^(x)-^- ) (y)-C JV )F(x ) y)=Q 
with the constant 

C N = - a 2 (l + k))N. 

It is important to note that also H^ N ' is a Schrodinger operator of the same 
kind as Hn, only the coupling parameters arc different. For the convenience of the 
reader we give the modified parameters for our examples I and III- VII, 

cij_i = — lu, c_i = — k — c, a_i = k — a, 6_i = k — b. 

The key property that makes the special case of Proposition ^. II stated in Corol- 
lary 12.21 interesting as a tool for constructing the eigenfunctions of a Schrodinger 
operator ([3|) is that the variables X of TL in (fT6|) are divided in two groups x and y 
such that the interaction terms involving variables belonging to different groups all 
vanish. It is interesting to note that there are other cases where this happens. For 
example, the number of particles with parameters mj = 1 and irij = — 1 need not 
be the same. More generally, if we divide the variables X into two groups where 
the 'mass parameters' in the first group are nij = 1 or —1/k and in the second 
group rrifc = —1 or 1/k, then all interaction terms between particles in different 
groups vanish (since \n(mj + mi,)(Kmjmk — l) = in all these cases). We then ob- 
tain an identity as in (|25p but now involving two deformed Calogero-Sutherland 
operators ^ which can have, in general, different particle numbers. Such an 
identity involves four different kinds of particles, and the number of particles of 
each kind can be arbitrary. These identities will allow us to construct eigenfunc- 
tions of the deformed Calogero-Sutherland operators ((6]) in Section |H The most 
general such identity corresponds to the following special case of Proposition 12.11 
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TV = N + N + M + M with non-negative integers N, N, M, M such that TV is 
non-zero, rrij — 1 and Xj = Xj for j — 1, 2, . . . , TV, tun+j = — 1/k and Xn+j = xj 
for J = 1, 2, . . . , N, m N+f j +k = -1 and X N+] ^ +k = y k for k = 1,2,..., M, and 

TO W+W+A/+A- = V« and ^jV+A>+M+K = WC for If = 1, 2, . . . , M. 

Corollary 2.3. Lef H N j^(x,i) = H N ^ be the operator in (|6|) urai/i V = V_i/ re , 

A/ M 

<~i(y>y) = E ( - C + ^*)) - £ «( - ^ + 

fc=l K=l 

+ «(« - !) E + E &) 

i<fc j<k 

+ (l- K )J2w( yj ,y K ) 
where M and M are arbitrary non-negative integers, and 

N N M M 

j = l J=l fc=l A=l 



Uj <k ( z k - z 3 r Y[ j<k {zk - zj) l,K 

X 



X 



U.j,j(Wj - Wj) 

U.j.k(™k - Zj) Y[j k {w k - ~zj) 

Ylj,k( w k - ZjY I\j,k( w k - z/) 1/K 

with Zj = z(xj), Wj — z(yj), zj — z(xj), wj — z(yj), and V m {x) and ipo,m(x) as 
in (|17| and (j!9|) . respectively, for m = — 1 and m = ±1/k. Then 

( H N,N^^ ~ H M?At&>y) - C N,N,M,Mj F N,N,MM(^^y?y) = 

with the constant 

Cmnm.m = ((N- - N_/ K ) 3 - TV_ + A>_ / K 3 



- -^-^ + ^) (( N _ ^/k) 2 -N + - N + /^ 
where 

N± = N±M, N±=N± M. 

Using this identity we can straightforwardly generalise our solution method 
and obtain many different formulae for eigenfunctions of the deformed Calogero- 
Sutherland operators © : for fixed particle numbers N and N one is free to choose 
M and M arbitrarily and each choice gives a family of eigenfunctions labelled by 
integer vectors n £ Z M+M . As discussed in Section g21 we obtain in this way a 
number of series representations for each eigenfunction, a fact which is interesting 
already in the standard (non-deformed) case N = 0. 
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2.3. Case II. In order to put the general discussion in the previous two sections 
into perspective we consider in some detail the special case II; see Tables [TJ and [2] 
The Schrodinger operator ([3]) is in this case given by 

N 



4 sin 2 



As we mentioned in the introduction, this Schrodinger operator was introduced 
and studied by Sutherland [50,51], and its eigenfunctions are given by the so-called 
Jack polynomials. In this special case the identity in Corollary 12.21 is equivalent 
to a well-known and important identity in the theory of Jack polynomials, and 
Corollary 12.31 gives a natural generalisation of this identity to the deformed case. 
More precisely, we have the following: 



Corollary 2.4. Let 



F(x,y) = [I 



. , 3 Uj,k( Wk ~ z ^ K 

with Zj = e lXj and Wj = e tyj . Then 

(26) (H N (x)-H N (y))F(x,y) = 0, 

where the arguments indicate that the Schrodinger operator acts in the variables x 
and y, respectively. 

As mentioned above, the identity (|2T)|) . and in particular the function F, is 
directly related with the theory of Jack polynomials. In order to make this precise 
we observe that 

F(z,w) = Y[(z k - Zj) K (l/ Wj - l/w k ) K tt(z : vs) 

j<k 



with the function II given by 

n(z,w) 



i 



It is well known that this latter function II has the following expansion in the monic 
Jack polynomials J\\ 

(27) n(z,w) = ^6 A J A (z)J A (w- 1 ) 

A 

for some (explicitly known) coefficients b\, and where w _1 = (1/wi, . . . , 1/iojv); see, 
e.g., Sections VI. 4 and VI. 10 in [39]. We recall that the monic Jack polynomials 
are usually denoted P\ rather than J\. We have used the latter notation in order 
to avoid clashes with notation used in later sections. Up to degeneracies in the 
spectrum of the Schrodinger operator Hn, this expansion can be deduced as a 
consequence of the identity ([26]). To establish this fact we let ^o(z) — Ylj <k (zk — 
Zj) K , and similarly for ^o( w1 )- It is well known that the functions ^>o(z)J'x(z) are 
eigenfunctions of H^; this can for example be inferred from results of Stanley [49]. 
Since -ff/v(y) is invariant under the substitution of —y,j for yj which maps Wj to 
1/iUj, the functions 5'o(w _1 )J7a(w _1 ) are eigenfunctions of the same Schrodinger 
operator Hn, and with the same eigenvalues. We observe that II(z,w) is invariant 
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under permutations of both the variables z as well as the variables w. It follows 
that it has an expansion in Jack polynomials of the form 

n(z,w) -^g A (z)j A (w- 1 ), 

A 

where the functions Q\ are to be determined. Since vE , o( w_1 )i7a(w _1 ) is an eigen- 
function of H^, the identity ([26]) implies that also \&o( z )Qa(z) is an eigenfunction of 
Hn, and with the same eigenvalue. It follows that, unless this eigenvalue is degen- 
erate, Qa(z) is proportional to J\{z). We mention that although such degeneracies 
do occur they are quite rare in this case. 

The identities obtained in the previous section thus suggest that there exist 
expansions, similar to the one given above, for each function F corresponding to 
a Schrodinger operator (|3|) or a deformed counterpart (|6|). However, in most cases 
the Schrodinger operator i?/v, as well as its deformed counterpart H N ^, is not 
invariant under a change of coordinates mapping Wj to l/wj. This explains the 
appearance of the operator and the constant Cat in Corollary 12.21 as well as 

that of H { M ] ~ { and C N fj M ^ { in Corollary [231 

3. ElGENFUNCTIONS OF CALOGERO-SUTHERLAND TYPE MODELS 

In this section we construct our series representations of the eigenfunctions of the 
Schrodinger operators ([3]) . These series representations are in terms of the functions 
/ n defined in ([7|). As we will see, the identities deduced in the previous section play 
a central role in this construction. These series representation generalises those 
obtained in [30] for the special case II and in [20] for I and IV (c.f. Table [1]). 

3.1. Construction of eigenfunctions. We fix the polynomials a, j3 in ([9]) and 

the coupling constant k > and consider the resulting Schrodinger operator Hn as 
defined by ([3]), ([TP]) . (fTTj) and ([13)) . We observe that the eigenfunctions in (f4]) are 
completely determined by the functions P n in N variables z, and we will refer to 
the latter as reduced eigenfunctions of H^. We also observe that P n is a reduced 
eigenfunction of if and only if it is an eigenfunction of the differential operator 

N , x 

(28) H N := ^(Hx-EoWo -jfafe^+WJ-fc^^, 

j=i &k z i Zk 

where the last equality follows by a straightforward computation using the fact 
that \?o is an eigenstate of Hn with eigenvalue Eq and the identities in (|HJ), ([T3|) . 
and (jllj) . Using the identity in Corollary 12. 21 it is now straightforward to compute 
the action of this differential operator on the polynomials / n . In order to state this 
action in a simple form, and also to facilitate a discussion of its implications, we first 
introduce some convenient notation. An important ingredient in our construction 
is the following partial ordering of integer vectors n, m E 7L N : 

(29) m ^ n O riij H h mjv < rij H ^n N , Vj = 1, 2, . . . , N. 

That this is only a partial ordering is easily seen, e.g., (522) and (441) are incompa- 
rable. The algebra of symmetric polynomials -P(z) in N variables z = {z\, . . . , z/v) 
with complex coefficients is denoted A at. The linear subspace consisting of the ho- 
mogeneous symmetric polynomials of degree n, together with the zero polynomial, 
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is denoted by Ajy. We also denote as ej the standard basis in Z N , i.e., (ej)k = Sjk 
for all j,k — 1, 2, . . . , N, and let 

E££ = (l-i/) ei + (l-p + i/)e fc . 

As before, 

(30) |n| := m + . . . + n N . 

For simplicity of notation we shall to each n £ Z N associate the shifted integer 
vector n + = (n^ , . . . , n^) with 

(31) n+ =n } +K(N + l-j). 
We are now ready to state and prove the following: 
Lemma 3.1. For each n £ *Z N , 

(32) H N f Il = (E I1 -E )f I1 

N 

- (am+(n+ - 1) + (A, - (1 + «)ai)(n/ - l))/ n -e, 
i=i 

JV 

2 oo 

+ k(k- l)^^^a p (2i/-p)/ n _ E P-. 
j<fc p=o i/=i 



Proof. Note that the function F in Corollary 1 2 . 21 can be written as follows: 

nv \ lTr / \ TT -«;(iv+i-j) n^fel 1 ~ w j/ w k) K 

F(x > y)=* (x)G(y)n^ fcfr _ Zj>fc)<c 

with 

G(y) :=n^o,-i(%). 

3=1 

The definition Q of the polynomials / n thus implies 

F(x,y)=* (x)G(y) ]T /„(z) W - n+ . 

It follows from Corollary 12.21 that 

(33) J] (i/ N / n (z))w- n+ = ]T /n(z)(^ + CW)w" n+ 

neZ" nGZ™ 

where 

:=G- 1 (y)^ ) G(y). 
Using ([8|), (fT3|) . and (fTT|) . as well as the definition of H N , we find that 

= -£ (aw^ + /m^,) + «(« - DE - A 

j=l j<k ^ ^ ■? ' ' 

with 

/3_iH = -p(w) + (1 + k)c/(u>). 
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Recalling ([9]) we expand the interaction terms in Laurent series and obtain 

a(wj) + a(u)k) 2a 2 WjWk + a\(wj + Wf.) + 2a 

a 2 = 



(Wj-Wk) 2 ~ W 2 (l - Wj/Wk) 



2 



wr 1 



p=0 v=l W k 

valid for \wn\ > \wn-i\ > ••• > \wi\ > maxk(\zk\)- It is now straightforward to 
compute the r.h.s. of (J33J), and by comparing the coefficients of w~ n+ on both sides 
of the resulting equation we obtain (|32[) with 

N 

E n = -^2[a2n+{n+ + 1) + (ft - 2(1 + ti)a 2 )n+} + C N . 

3=1 

By straightforward computations we find that the latter coincides with E n in QI4[) 
and SB. □ 



It is clear that the action of Hjq on the polynomials / n is triangular in the 
sense that Hn / n is a linear combination of polynomials f m with m ^< n. In fact, 
depending on the specific choice of a and /3, the set of possible m is in many cases 
significantly smaller. In any case, this triangular structure suggests that there exist 
eigenfunctions, with corresponding eigenvalues E^ — Eq, of the differential operator 
Hn which are of the form 

(34) P„ = /„ + E«n(m)/m 

m 

where the sum is over integer vectors m -< n. Indeed, if we take this as an ansatz 
then we find that this is the case if and only if the coefficients satisfy the recursion 
relation 

N 

(35) (E n - E m )u n (m) = - ^2m+(f3 + a\(mj - /c))it n (m + e,-) 

3=1 

N 2 

- a J^m+(m+ + I)u n (m + 2e^) + k(k - 1) },/^ a p (2v - p)u n (m + E^). 

j=l j<kp=0 

Since u n (n) is fixed to one, it is clear that if E n — E m ^ for all m which appear 
in (|34p , then this recursion relation uniquely determines the remaining coefficients 
it n (m). This condition of non-degeneracy on the eigenvalues is generically satisfied, 
but there exist special cases and specific integer vectors n for which it fails to hold 
true; see Remark 13.21 for a further discussion of this point. 



3.2. Case II. In order to draw attention to some of the specific features of the 
action of Hn on the polynomials /„ we consider in some detail the special case II; 
see Tables Q] and [2j On the one hand this is in many ways the simplest non trivial 
case, and our construction of reduced eigenfunctions is therefore particularly simple 
in this case, but, on the other hand, it already contains many of the key features of 
the general case. For the special case II we have a(z) = —z 2 and (3(z) = — z, and 
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consequently that 

oo 

H N f n = (E n - E )f n - 2k(k ~ 1) £ £ v fa-v(e h -e } )- 

j<kv=l 

Clearly, the right hand side now only contains polynomials / m such that |m| = |n|. 
It is interesting to compare this expression with the corresponding action of Hn on 
the monomial symmetric polynomials 



(36) m x (z):=J24 ni) 



•VPCI) „ X P(N) 

Z N 



with A a partition of length t(X) < N, and where the sum extends over all distinct 
permutations P of the parts Xj of A. To this end we let [a; J denote the integer part 
of x G K. In addition, for any n £ Nq we identify m n with the unique ui\ such 
that n = -P(A) for some permutation P. Using the identity 

n—m 

— (x 2 d x - y 2 d y )(x n y m + x m y n ) = (n - m) E x n - k y m+k + 2mx n y m , 



x-y 



k=0 



valid for any non-negative integers n and m such that n > m, it is now easy to 
verify that 

H N m\ = (Ex - E )m\ + 2k 2j(Aj - A&) 2J m \+v(e k - ej )- 

j<k v=Q 

We note, in particular, that the action of Hn on the polynomials / n does not 
contain an explicit dependence on the integer vector n. In contrast, the action on 
the symmetric monomials m\ do contain an explicit dependence on the partition 
A. It is also interesting to note that the corresponding recursion relation ([55)) 
can in fact be solved to yield the following explicit series representation for the 
eigenfunctions of H^: 



^„ = /„+£(2K(i-K)r Yl E^ x -" x E 

s=l ii<fei! y i = l js<k s u s = l 

s _j 

x II ( Ex ~ E ^-J2UrM^ t -B h )) fx-T, s r =iM^ r - ejr h 

r=l 

see [31] for further details. It is important to note that the fact that /„ is non-zero 
only if n ^ (see Corollarv l3.ip implies that this series only contains a finite number 
of terms, and thus is well-defined. Even though it can be done in principle, it 
seems harder to solve the recursion relations resulting from the use of the monomial 
symmetric polynomials. Detailed discussions of the above construction of series 
representations for the reduced eigenfunctions of the Schrodinger operator ([3]) in 
the cases I and IV can be found in [17,20]. We also mention that it is possible 
to write down explicit formulae for the reduced eigenfunctions of the Schrodinger 
operator ([3]) in the general case which, however, are somewhat involved^ 



The interested reader can find these formulae in the first arXiv version of the present paper; 
see http:/ /arxiv.org/abs/math-ph/0703090vl 
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3.3. Completeness of the reduced eigenfunctions. In Section 13.11 we con- 
structed, under a certain condition of non-degeneracy on the corresponding eigen- 
values, for each integer vector n £ Z N a reduced eigenfuntion P n of the Schrodinger 
operator H^. In this section we show that if we restrict attention to integer vectors 
n = A for some partition A, then the corresponding reduced eigenfunctions form a 
linear basis for A at. The idea of the proof is to first establish that the expansion of 
the symmetric polynomials / n in terms of a particular set of symmetric polynomials 
g\, known to constitute such a basis, has a particular triangular structure. This 
implies that the f\, labelled by the partitions A of length ^(A) < N, form a linear 
basis for Ajv- We then obtain a similar result for the reduced eigenfunctions P\ by 
showing that they are, in the same sense, triangular in the fx- 

We proceed to recall that the polynomials g\ in question, sometimes referred to 
as the modified complete symmetric polynomials, can be defined by the expansion 

(W- 1 = Of 1 ,...,!^ 1 )) 

(37) n(z,w) =^. 9A (z)m A (w- 1 ) 

A 



valid for minfc(|i0fe|) > maxj ( | Zj | ) , and where the sum extends over all partitions 
A of length £(A) < N. It is well known that the g\ constitute a linear basis for 
An; see, e.g., Section VI. 10 in [39]. Using this fact we now establish the following 
relation between the / n and the g\, which, in particular, implies that the / n indeed 
are well-defined symmetric polynomials: 



Lemma 3.2. Let n e Z N . Then 



fn — / a npiff/x 



for some coefficients a n /x and where the sum is over partitions fi -< n such that 
\fjt\ = |n|. Moreover, if n = A for some partition A of length £(X) < N, then 
a\\ = 1. 

Proof. It follows from the definitions of the g\ and / n that 
£ / n (z)w-" = [](l-^M-) K E 5 P (m)(z)w- m 

OO / s 

= e 5p ( m)(z)n e^ 1 ) 9 ^ K w- m+ ^<*«*^-*) 



mGN™ j<kq jk =0 



where p(m) denotes the unique partition obtained by permuting the components 
of the integer vector m, e.g., p(2,4, 5,2) = (5,4,2,2). Comparing the coefficients 
of w~ n on both sides of this identity we find that 

(38) /„ = [] E (-V* h L. t Up(n +Ei< ,^(e i -e fc )). 

We thus obtained a representation of / n as a linear superposition of polynomials 
g^ where fj, = p(m) for 

m = n + £ qjkiej - e fe ) 

j<k 

with non-negative integers q^. Obviously, |ju| = |m| = |n| for such fi. To complete 
the proof we thus only need to show that /i^n. For that we observe that m ■< n. 
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Moreover, since /x is a partition, ji\ > fj,2 > • • • > 1-J.n, and the definition of the 
ordering (|29|) thus implies that /j, ^ Pfi for any permutation P of the parts of fi 



Remark 3.1. It is interesting to note that in the particular case k = 1 the result 
in (|3"5)) is equivalent to the so-called Jacobi-Trudi identity (originally due to Jacobi 
[23]), which implies that the well-known Schur polynomials s\ are identical with 
the polynomials fx for k = 1. The proof of this fact can be obtained as a simple 
consequence of well-known arguments which can be found in Macdonald's book [39] ; 
see the proof of (3.4") and the preceding discussion in Chapter I. 

The fact that each g^ is a symmetric homogeneous polynomial of degree now 
implies the following: 

Corollary 3.1. The function f n , n € Z N , is non-zero only if n y 0. In that case 
it is a symmetric homogeneous polynomial of degree |n|. 

Let n be a non-negative integer and consider the two sets of polynomials g\ and 
fx which are both indexed by those partitions A of n which are of length £(X) < N. 
Let K = (Kxy) denote the transition matrix from the fx to the gx, defined by the 
equalities fx = Kx^g^- Given a partition A, Lemma l3~2l implies that Kx^ = 
unless ^ A and that Kxx = 1- L is readily verified that the inverse of such a 
matrix always exist and is of the same form; see, e.g., Section 1.6 in [39]. It follows 
that 



where (K 1 )xfj, are the elements in the matrix inverse to K and the sum is over 
partitions fj, -< A of n. We thus obtain the following: 

Proposition 3.1. For any non-negative integer n, the polynomials fx, labelled by 
those partitions A of n which are of length £(\) < N, form a linear basis for A^-. 

Suppose now that the reduced eigenfunctions P n of the form (f34|) exist for all 
integer vectors n S 2j N such that n = A for some partition A. From Lemma 13.21 
and (|39p we have that any such reduced eigenfunction has a series expansion of the 
form 



for some coefficients vx^, and where the sum is over partitions jU -< A. Fix a 
partition A of length ^(A) < N and consider the transition matrix from the reduced 
eigenfunctions P^ to the polynomials f ^ , both indexed by the partitions fi of length 
^(m) ^ N such that p ^ A. Observe that its entries are non-zero only if fi ^< A and 
that all diagonal entries are equal to one. It is clear from the arguments leading 
up to Proposition 13.11 that this transition matrix has a well-defined inverse of the 
same form. It follows that each polynomial fx is a linear combination of reduced 
eigenfunctions P^ with fi ^ A. Proposition 13.11 thus implies the following: 

Proposition 3.2. // the reduced eigenfunctions Px in (|34fl exist for all partitions 
A of length £(X) < N then they form a linear basis for A^r. 

Remark 3.2. As discussed after Equation (|35|) . a reduced eigenfunction of the form 
(|34| exist if E n — E m ^ for all m which appear in (|34|) . It is easy to see that 



and, in particular, fi ^ m. 



□ 



(39) 
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this condition of non-degeneracy holds true if a 2 = (and (3\ ^ 0), including the 
special cases I and IV (see, e.g., [20]). It is also known to hold for all integer vectors 
n G Z N such that n = A for some partition A in the special case II (see, e.g., [31]), 
assuming that k > 0. It is interesting to note that if a 2 7^ and if one allow k 
to be negative then there exist parameter values and integer vectors n such that 
E n = E m for some integer vector m which appear in (|34|) . We stress, however, that 
for generic parameter values and integer vectors n this condition of non-degeneracy 
on the eigenvalues is satisfied. 



4. ElGENFUNCTIONS OF DEFORMED CALOGERO-SUTHERLAND TYPE MODELS 

In this section we sketch the generalisation of our construction of eigenfunctions 
from the previous section to the deformed Calogero-Sutherland operators H N ft (p]). 
emphasising the features which are new. Just as the construction in the previous 
section was based on Corollary [221 the discussion below is based on Corollary 12. 31 
We recall that the identity in this latter corollary relates the operator H N ft to an 
operator of the same type, but with particle numbers M and M. We stress that 
these particle numbers can be chosen freely, and that, as we will discuss in Section 
14.21 the complexity of the resulting series representation of a given eigenfunction can 
be highly dependent on the specific choice. As described already in the beginning 
of Section [TJ this can be used to construct simpler explicit formulae for many of 
the eigenfunctions already in the special case N = 0. 

We proceed to describe the structure of the eigenfunctions of the deformed 
Calogero-Sutherland operators Following Sergeev and Veselov [48] we let 

A N ft K be the algebra of polynomials P(z, z) in two sets of independent variables 
z = (zi, . . . , zn) and z = {zi, . . . , Zft) with complex coefficients such that they are 
separately symmetric in the variables z and z and, furthermore, obey the condition 

{d Zj +Kd Sj )P(m,x)\, i=Sj =0 

for all j — 1, . . . , N and J = 1, . . . , N. We will show that the deformed Calogero- 
Sutherland operators ^ have eigenfunctions which are naturally labelled by par- 
titions A = (Ai, A 2 , . . .) in the so-called fat (N, iV)-hook, i.e., Xjy+i < N, and are 
of the form 

(40) *a(x,x) - * (x,x) J P A (z,z) 
with Zj — z(xj), zj — z{xj), 

(41) *o(x,x) = [[M x j> 11 ipa,-i/ K (xj) ff — j-~ ^ , 

3=1 .1 1 llj,A Zj Z J> 

and where the P\ are particular polynomials contained in the algebra ft . 
Moreover, the corresponding eigenvalues are given by 

(42) E X = E -J2 (aiXjiXj - 1) + (2a 2 ( K (N - j) - N) + ft) Xjj 
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with 

(43) E ^-^-^^N ~N/k) 3 -N + N/k 3 

- + ( - N/n) 2 -N- N/^ . 

To the best of our knowledge, this fact was previously known only in special cases; 
see [47,48]. 

4.1. Construction of eigenfunctions. At this point we fix the polynomials a 
and f3 in and the coupling constant k > and consider the resulting deformed 
Calogero-Sutherland operator H N jy © . As we will see, the construction of explicit 
series representations for the polynomials P\ is very similar to the one in the pre- 
vious section, but there are a few important differences in definitions and notation 
which we now specify. 

We start by noting that ty\ is an eigenfunction of H N $ if and only if P\ is an 
eigenfunction of the differential operator 

N 

(44) H N fj := % \H N j - E )* = - £ («(^ )^ + P(z 3 )d Z] ) 

+ k £ (a(z J )9| J + ^ 1/K (zj)d- Zj ) - 2 K ^2^-d Zj 
j=i &k z i Zk 

+ 2 y «(^) g _ 2K y _2&L dii - 2 V a ( g *>a,. 

with corresponding eigenvalue — £-o! this equivalence follows from a straightfor- 
ward computations using that (HO is an eigenstate of H N ^ © with eigenvalue 

E (Corollary [231 for M = M = 0) and the identities ©, (13]), and |TT]). We shall 
refer to the polynomials P\ as reduced eigenfunctions of the deformed Calogero- 
Sutherland operators ©■ Furthermore, we let 



(45) F^ M ^(z,z,w,w) 
T lj<k 



n i<fc (i - "iM)"i]j< K (i-"j/%) 1/B 



rij.fcCi - z j/ w k) K rwa - 

where 

(z,z) = (zi,...,2jv,2i,...,%) and (w, w) = (wi, . . . , wat, u>i, . . . , w^) 

and define the polynomials f^ ,M \z,i), n 6 Z M+A/ , through the following expan- 
sion: 

(46) F w ^ MJ ^(z,z,w,w) 



™i . . . w ~ nM iTr nM+1 

W M w l ~ M 
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valid in the region \wi\ > ■ ■ ■ > \wm\ > \wt\ > ■ ■ ■ > \w^ f \ > m&Xj.j(\zj\, \zj\). The 
next step is to compute the action of the differential operator H N ^ on these poly- 
nomials f^ [ ' M \ In order to simplify such a computation it is useful to introduce 
the following 'parity' function q on the index set {1, . . . , M + M}: 

(47) q(j) := 



0, ifj 1 \/. 

1, i£j = M + l,...,M. 



We also find it convenient to define the following shifted integer vectors n + associ- 
ated with each quantum number n: 



(48) n 



nj+K(N+l-j)-N, if j = 1,...,M, 

rij + (N + M + 1 - j) /k-N + M, if j = M + 1, . . . , M. 



Proceeding in analogy with the proof of Lemma 13.11 in the previous section it is 
now straightforward to obtain the following: 

Lemma 4.1. For each n e Z M+M . 
(49) H^ffM = (2#*.*> - Eo) 

M+M 

((-^K - 1) + (A - (l - i-K) 1 -*®) ai) (nf - 1)) 



3 = 1 

M+M 

r(M,M) 



-ao E(--) 9(i) (<-l)K+-2)/^ 



2 oo 

9(j)-9(*) V" \^ n _ (0,, - n \ f( M <™\ 



j'<fe p=0 i/=l 



-e; 



with 



M 



(50) £^ M > M ) - So = - J2 («2^(^ - 1) + (2a 2 ( K (N - j) — N) + fa) nj ) 



3 = 1 

M+M 



+ J2 (na2n J (n j + l) + (2a 2 {N~M + l-j-K(N-M))-f3 1 )n J ). 

j=M+l 

Similarly to the previous section, Lemma 14.11 implies that that there exist an 
eigenfunction of the form 

(5i) f* u >*> = f^+j2M^)f^\ 
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where the sum is over integer vectors m £ jm+M sucn ^hat m ^! n, if and only if 
the coefficients u n (m) satisfy the recursion relation 



J £ (M,M)_ £ (A/,M)j Hn(m) 



M+M 



£ m+ ((-/e)«0>ai(n+ + 1) + (fa - (l - (-k) 1 -^) ai)) Un (m + ej ) 

M+M 

- a (-K) g(j) mj(mj + l)u„(m + 2ej) 



2 OO 



+ (1 - «) £(-«)i-<Ki>-*W £ £ ap(2l/ _ pK(n + E ^). 

j'</c p=0 l/=l 

Provided that E n M ' M ' — E^ ^ for all integer vectors m which appear in (151]) 
this uniquely determines all coefficients u n (m). As previously mentioned, an im- 
portant feature of our construction is that the non-negative integers M and M can 
be chosen freely, and for each choice we obtain a family of reduced eigenfunctions 
p(M,M) q £ jj which are labelled by integer vectors n 6 Z M+Af . This yields 

eigenfunctions (|40|) of H N ^ with corresponding eigenvalues En given by (|50)) 
and (|4"3")) . A priori, it is not clear whether these eigenvalues coincide with the ones 
given in (|42[) . In particular, it seems that the eigenvalues we obtain depend on 
the specific values of M and M, but, of course, the spectrum of H N ^ does not. 
This apparent paradox can be resolved by exhibiting an explicit mapping from the 

eigenvalues E n M ' M ^ (|50|) to E\ (|42[) under certain restrictions on the integer vec- 
tors n G Z M+M . To state this result we need the notion of the conjugate fj,' of 
a partition /x, obtained by interchanging rows and columns in the Young diagram 
corresponding to /x, e.g., 



M = (5,3,2) 



see, e.g., Section I.I in 



(3,3,2,1,1); 



Lemma 4.2. Let n = (m, fi) with m £ Z M and fi a partition of length £(fi) < M . 
Then 

(52) = E x 

with A = (m, fj,'). 

Proof. We recall the following two well-known identities obeyed by all partitions \i: 



where the sums are over all non-zero parts fij and fj,'j of fi and fi', respectively. 
This, together with the identity ^ ■ Hj = J2j inserted into the second sum in 
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([50]) . yields 



M 

E 



(a 2 m j (m j - 1) + {2a 2 {n{N - j) - N) + (3i)mj) 



£( M ') 

E 



{a 2 n'M - !) + ( 2 M><^ - M - - N) + ft)^). 



The statement now follows by substituting Aj for and Am+j for /i'-. 



□ 



4.2. On the complexity of the series representations. The construction de- 
scribed in the previous section gives reduced eigenfunctions P n for all integer vectors 
n such that a certain condition of non-degenracy on the corresponding eigenvalues 
is satisfied, but unless n satisfies the condition in (|52p for some partition A in the 
fat (N, 7V)-hook this function ought to vanish identically. Lemma l4~2l gives . for each 
partition A in the fat (N, iV)-hook, such an integer vector n provided one chooses 
(M, M) = (N,N), and the latter choice thus is natural if one is interested in a 
general formula which works for (essentially) all eigenfunctions. However, we stress 
that we have constructed a family of reduced eigenfunctions for each choice of non- 
negative integers M and M , and that there exist a simple relation between any two 
such families: For any two sets of non-negative integers (Mi, Mi) and (M 2 ,M 2 ) 
one can show that two reduced eigenfunctions pU^i, M i) anc j p^' /2 ' M2 ^ labelled 
by integer vectors n! and n 2 which correspond to one and the same partition A 
under the mapping defined in Lemma 14.21 are equal up to normalisation; see the 
first arXiv version (vl) of [18]. We can thus, in a simple manner, obtain different 
series representations for one and the same reduced eigenfunction by varying the 
values of M and M. Since the complexity of this series representation is highly 
dependent on the values of M and M, one can choose the latter such that the com- 
plexity is minimised, in many cases substantially below that of the canonical choice 
(M, M) = (N, N) . This is already evident in the case of the Schrodinger operators 
([3]), i.e., for N = 0. For example, suppose that we are interested in the reduced 
eigenfunction of the standard Sutherland model for N = 8 and corresponding to 
the partition with diagram 



Rather than setting M = 8 and M — we can use the fact that the conjugate of 
this partition is given by the diagram 



and thus the sought after reduced eigenfunction can equally well be obtained by 
setting M = and M = 2, in the process decreasing the complexity of its series 
representation from the 8-particle case to that of only 2-particles. As a further 
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illustrative example we consider the reduced cigenfunction corresponding to the 
partition with diagram 



for which we can set M = M = 2 to minimise the complexity of its series represen- 
tation. Similar examples can be easily constructed also in the 'deformed' case. In 
generally, it can be readily verified that, given a specific partition A, the minimal 
value for M + M is attained by setting M — j and M = Xj+i where j denotes the 
row in the diagram of A for which j + Xj+i is minimal. This means, in particular, 
that we always can choose M and M such that their sum does not exceed the length 
£(A) of the partition A in question, which may be less than N + N. This observation 
is also reflected in the fact that fn M ' M ^ = f^' MJrK > f or a \\ non-negative values of 
K as long asM + M > £(n); the latter is readily inferred from the definition of the 
polynomials /n'^. 

4.3. Elementary properties of the reduced eigenfunctions. We conclude this 
section by verifying that the functions fn*' M \ and consequently also the reduced 
eigenfunctions Pn' M \ are contained in the algebra Ajy ^ for all values of the 
non-negative integers M and M. 

Lemma 4.3. Let n e Z A/+ * / . Then the function fn is non-zero only ifn >z 0, 
and in that case it is a homogeneous polynomial of degree |n| in A N ^ . 

Proof. The proof is by straightforward computations: expanding each factor on the 
r.h.s. in (|45[) in the region \wi\ > ■ ■ ■ > \w M+ ^ I \ > max^jQ-Zjl, \zj\) (recall that 
Wj = wm+j) in a binomial series one obtains a well-defined series representation of 
F N fj M £j as a linear superposition of monomials w~ n . By inspection one can check 
that in this series only terms w _n with n y appear, that each such term appear 
only a finite number of times, and that the coefficients of w~ n is a homogeneous 
polynomials of degree |n| in the variables Zj and zj. It thus follows from ([41))) that 
j(m,m) ^ non _ zer0 only if n y 0, and in that case it is a homogeneous polynomial 
of degree |n|. To prove that the /n M ' M ' ) are contained in A N we observe that 
F N M £ { is separately symmetric in the variables Zj and zj, and that 

(d Z] + ^ ij )F NfjM ^ I \ z _ ij = 

for all j = I, . . . , N and J = I , . . . , N. The latter follows from (|45|) by straightfor- 
ward computations, 

M j M - 

a Zj logf^ v u u = — - — k ^ 
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and 



9zj log Fn,n,m,m 



e- 5 — -i:^-^- 

L — ' Z J — Wh K * — ' Z J — Wk 

fc=l ■' K K=l J " 



□ 



5. Concluding remarks 

In this final section we briefly discuss the relation of our construction to the 
so called generalised classical orthogonal polynomials of Lassalle and Macdonald, 
present some further remarks of the effect of spectral degeneracies on our construc- 
tions of reduced eigenfunctions of the Schrodinger operators Q and the deformed 
Calogero-Sutherland operators ©, and comment on the question of completeness 
in the case of the deformed Calogero-Sutherland operators. 

5.1. The relation to generalised classical orthogonal polynomials. As is 

well known, the spectra of the Schrodinger operators ([3]) are degenerate. In partic- 
ular, in cases I and IV the eigenvalue of a given reduced eigenfunction P n depends 
only on the weight |n| of the integer vector n. An interesting question is thus 
what the relation is between the basises for the reduced eigenfunctions we con- 
struct and previously known such basises. Below we will sketch a proof of the 
fact that our basis coincide, up to normalisation, with the generalised classical or- 
thogonal polynomials of Lassalle [33-35] and Macdonald [38] . We recall that these 
latter polynomials are reduced eigenfunctions of the Schrodinger operator in case I 
(Hermite), case IV (Laguerre) or case VI (Jacobi), and of the form 

for some coefficients ux^, and where /jCA means that the Young diagram of fi is 
contained in that of A. We shall require the following partial order on the set of 
partitions: 

H < A m -\ 1- fXj < AH h Xj, Vj. 

When restricted to the set of partitions of a fixed weight this coincides with the 
so called dominance order. It is a simple exercise to verify that if \fi\ = |A| and 
(j, < A then fi >z A and vice versa. Furthermore, comparing the expansions (|27jl and 
(|37| , and using the triangular structure of the Jack polynomials when expanded in 
monomial symmetric polynomials, one finds that 

.9a = ^ aXfiJ/j, 

— | A| ,/x> A 

for some coefficients such that a\\ = b\; see Equation 10.16 in Section VI. 10 
of [39] for the definition of b\ (in doing so observe that the parameter a = 1/k). 
It is now readily inferred from Lemma 13.21 that any reduced eigenfunction P\ , as 
constructed in Section l3~H is of the form 

/x^A 

for some coefficients vx^ such that v\\ = b\. Since /i ^ A if /j C A and, up to 
normalisation, there can be only one eigenfunction of this form, it follows that the 
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relation between the reduced eigenfunctions we construct and those of Lassalle and 
Macdonald is given by 

P\ = bxVx- 

5.2. Spectral degeneracies. We stress once more the issue of degeneracies of the 
eigenvalues of the Schrodinger operators ([3]) and the deformed Calogero-Sutherland 
operators ([6]): since our construction of reduced eigenfunction of these operators 
only is valid under a certain non-degeneracy condition on the corresponding eigen- 
values, this issue plays a decisive role in our construction. From our point of view, 
it is essentially only this issue which truly distinguishes the different special cases 
of our results listed in Table [1] It thus would be interesting to investigate this issue 
further. 

5.3. Completeness for deformed Calogero-Sutherland operators. We men- 
tion that one can prove an analogue of Proposition 13.21 for our construction of re- 
duced eigenfunctions of the deformed Calogero-Sutherland operators ((6]). In fact, 
proceeding similarly to Section 13.31 one can prove that the reduced eigenfunctions 
constructed in Section [4.11 have a triangular expansion in so-called super Jack poly- 
nomials. We recall that the super Jack polynomials are known to constitute a 
linear basis for the algebra A N ^ . Using these two facts above one can show that 
if M > N and M > N then the reduced eigenfunctions constructed in Section |4~T1 
are also a linear basis for A N ^ K . If either M < N or M < N then they span only a 
subspace of A N ^ which, however, can be given a rather simple characterisation. 
A detailed discussion of these facts can be found in the first arXiv version (vl) 
of [18]. 
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